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Abstract 

A single molecule break junction device serves as a tunable model system for probing the many 
body Kondo state. The low-energy properties of this state are commonly described in terms 
of a Kondo model, where the response of the system to different perturbations is characterized 
by a single emergent energy scale, ksTx- Comparisons between different experimental systems 
have shown issues with numerical consistency. With a new constrained analysis examining the 
dependence of conductance on temperature, bias, and magnetic field simultaneously, we show that 
these deviations can be resolved by properly accounting for background, non-Kondo contributions 
to the conductance that are often neglected. We clearly demonstrate the importance of these non- 
Kondo conduction channels by examining transport in devices with total conductances exceeding 
the theoretical maximum due to Kondo- assisted tunneling alone. 
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Strong interactions between charge carriers can lead to the emergence of novel many-body 
transport phenomena. In some cases the resulting states may exhibit universal characteris- 
tics that bridge seemingly diverse systems. The Kondo effect is an archetype of such corre- 
lated electron behavior, in which a collective state is established when itinerant conduction 
electrons antiferromagnetically screen a local magnetic moment. 

Innovative device configurations possessing a small island tunnel-coupled to a Fermi sea 
allow for the formation of a localized magnetic impurity, represented here as a single un- 
paired spin. Examples include semiconductor quantum dots,|lL M nanowires,y] carbon 
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nanotubes,|4], |5], individual molecules, p, |7J including C 60 ,|8HlfJ and in STM measurements 
of both magnetic adatoms on metallic surfaces ll| and in chemically homogenous magnetic 
materials [l2j. Kondo screening occurs when the unpaired spin forms a singlet state with 
the conduction electrons as they tunnel on and off the island via elastic cotunneling pro- 
cesses. This results in the quenching of the local spin, and produces an extra peak in the 
density of states at the Fermi level (Figs. UK and b). The hallmark signature of the fully 
screened spin-| Kondo effect in these structures is an enhanced zero-bias conductance that 
forms below an emergent energy scale known as the Kondo temperature, Tk (&b = 1), which 
represents the binding energy of the many-body singlet state. 

The spin-i Kondo model describes the antiferromagnetic interaction between the mag- 
netic moment and conduction electrons at strong coupling in a particle-hole symmetric 
device jsl, [k], [jjj]. A remarkable feature of this model is the formation of a local Fermi 
liquid state characterized by a single energy scale, Tk- As a result, the low-energy behav- 
ior is expected to be universal. Realistic devices can deviate from this expectation for a 



number of reasons. This could occur, for example, due to a departure from particle-hole 



symmetry 



151 ] or due to multiple energy levels on the island participating in transport. To 



fully characterize the electronic properties of devices, a reliable extraction of Tk, and ad- 
ditional parameters, is critical to the ensuing analysis. Previous analyses for Tk and a 
number of transport coefficients, based on semi-phenomenological formulas for the spin-^ 
Kondo model, have found agreement between different gated semiconductor systems, but 
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systematic disagreement with single molecule devices. 

In this paper we report an analysis of conduction measurements through single molecule 
devices in the low energy regime as a function of temperature, bias, and magnetic field. 
We observe zero-bias resonance peaks with total conductance, G, exceeding the unitary 
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limit (Go > 2e 2 //i), which serves to highlight the presence of extra transport channels 
in addition to Kondo scattering. Our analysis allows us to characterize the low-energy 
behavior of these devices in a consistent manner in terms of an Anderson model and a 
background contribution. Most investigations ignore non-Kondo conductance channels 
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181 ]. as they are assumed to be vanishingly small or simply irrelevant. We demonstrate that 



this background contribution can be non-negligible in single molecule devices fabricated 
with a break junction technique. Properly separating Kondo and non-Kondo transport is 
essential in obtaining consistent values for Tk and the corresponding transport coefficients, 
and may impact our understanding of the effective low-energy model of real devices. 

Fits to the frequently employed phenomenological equation for G(T), together with re- 
cently developed analytical descriptions for G(V) and G(B) in the spin-- Kondo model (see 
below), provide three independent means of inferring the characteristic energy scale Tk- We 
determine a background conductance offset by including it as a free fitting parameter for the 
three fits, assuming that it is weakly dependent upon temperature, bias voltage and mag- 
netic fields for energies < ^bTk and may thus be modeled as a constant. The "optimum" 
offset due to background conduction is the one that most nearly produces the same Tk for 
the different perturbations. This treatment is applied to both new data and to previously 
published results, potentially resolving the conflicting coefficient values found in different 
systems for the ostensibly universal transport coefficients. 

The metallic break junction devices we used were defined by e-beam lithography on n + 
Si substrate with 200 nm SiC>2. They were composed of 15 nm Pd, deposited by ebeam 
evaporation. The resulting constriction patterns possessed minimum widths between 40 nm 
and 80 nm. Devices were cleaned first with a standard solvent rinse followed by an oxygen 
plasma to remove residual organic compounds. A single drop of 100 /iM Cgo in toluene 
solution was placed on an array of devices for 30 sec and then blown off with N 2 gas. Mea- 
surements were performed inside of a 3 He- 4 He dilution refrigerator with a base temperature 
< 100 mK. Thermally-assisted electro migration was employed at temperatures below 2 K 
to create a tunneling gap.fl9| Differential conductance (dl/dV) was measured as a function 
of source-drain bias (V) using standard low frequency lock-in techniques (Fig. Hb). Identifi- 
cation of devices possessing the desired electrode-C 60 -electrode geometry is discerned from 
the resulting transport characteristics, as has been discussed elsewhere |S|, |20 | 

The Kondo ground state is expected to be a Fermi liquid, and thus the lowest order 
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correction to G(T = 0, V = 0, B = 0) is quadratic in temperature, bias voltage and magnetic 
field. When normalized by Tk, the response of the system to different perturbations can be 
described by equations of the same form at low energies (kgT, eV, \q\hbB <C k B T K ), 



G(T) 


= G(T,0,0) 


= G(0,0,0) 


-c T G {T/T K ) 2 , 


(1) 


G(V) 


= G(o,v,o) 


= (7(0,0,0) 


- c v G (eV sd /k B T K ) 2 , 


(2) 


G{B) 


= G(0,0,B) 


= G(0,0,0) 


- c B G (\g\fi B B /k B T K ) 2 , 


(3) 



where G = G(0, 0, 0) — Gb with Gb being the background contribution and ct(c v , c b ) 
is a transport coefficient in T (V, B), respectively. The equilibrium conductance in the 
Kondo regime measured in a variety of systems has been shown to obey a universal 

temperature dependence described by the following phenomenological expression 18| that 



smoothly connects the low (T/Tk <C 1) and intermediate energy (T/Tk ~ 1) regimes, 



G(T) = G 



(2--!)(0 



Gh 



(4) 



where s = 0.21 for a spin-| Kondo model in the absence of any potential scatterer. [18|, |21| 
The Kondo temperature is defined by [G(T = Tk) — Gb]/Go = 1/2. Eqn. (BJ has been 
widely used as the principle means of estimating Tk in experimental systems. However Gb 
is often not included. 

We have measured several samples that exhibit Kondo conductance with total G above 
the unitary limit of 2e 2 /h. The zero-bias peak in sample A (FigsHH,e and|2^L,b), for example, 
exhibits typical Kondo resonance behavior. However, the total conductance cannot result 
entirely from cotunneling via the same single particle level. This implication of Gb > is not 
surprising in the case of break junctions, where the interelectrode gap is frequently < 1 nm, 
and the tunneling current is exponentially sensitive to the electrode-molecule displacement. 
Besides parallel conduction due to direct tunneling between electrodes, some non-resonant 
source of conduction is the postulated origin for the Fano-lineshape of Kondo resonance 
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peaks found in some break junction [12j and STM experiments. 

It has been proposed that data may be treated without accounting for extra conduction 
channels, la, |23| This can lead to coefficients closer to theoretically proposed values, but 
it also results in poor fits to data (Fig. top) and erroneous estimations of Tk for high 
conductance samples inconsistent with the commonly used approximation that the resonance 
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FIG. 1: (color online), (a) Energy level diagram showing charging energy, U, width of single 
particle energy level, T, depth of energy level, eo, below the equilibrium Fermi level, fj,, and the 
co-tunneling exchange process leading to Kondo resonance resulting from (b) an extra peak in the 
density of states at the equilibrium Fermi level, (c) Schematic of measurement setup, (d) G(V) 
vs. V for sample A at T/Tk < 1-5. (e) dl/dV vs. T also for sample A. Solid black lines are fit to 
Eqn. with Gf, = (top), Gb equal to optimum value as determined in text (middle), and Gb 



determined through procedure used in Ref 



171 ] (bottom). Gray dashed line is a fit to the analogous 



equation for an underscreened spin-1 Kondo effect (See Supplementary Information). 



peak's full width at half maximum is roug hly 2y/2k B T K [22|. In Ref. fl7j similar devices 
were studied and a background contribution to the overall conductance, linear in applied 
bias, was subtracted from conductance measurements. It was established via a fit to the 
trace dl/dV vs. V, using a range of applied bias, ±V, which incorporated both minima 
adjacent to the zero-bias resonance. As an example, this "full" background subtraction is 
shown at the bottom of Fig. [4^ when applied to sample A. 

Rather than relying solely on an estimate for Tk from Eqn. (J4]), we additionally utilize 
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expressions for G 



parameter. 
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V) and G(B) while using a constant conductance offset, Gb, as a fitting 
24j In Ref. 25J an expression for the zero temperature conductance out of 
equilibrium, G(V, 0, 0) for the spin-| Kondo model was proposed. This analytical form 
provides a description of the Kondo resonance peak analogous to Eqn. (J4j) in which 

(2 1 ^-l)x 2 



G(V) = G G 



«2 



(5) 



l + b(\x 

where x = V/T* K and (T K /T K ) 2 « tt. 

The behavior of the equilibrium conductance as a function of magnetic field can also be 



used to characterize the device. Extending the method of Ref. [15| to finite magnetic field 
yields a systematic expansion of the conductance G(T, V, B). For the spin-| Kondo model, 
one finds 

G(B) ~ ^ (6) 



1 ' V Ak B T K J 

where £ = 3T L T R /(T R + T L ) 2 is a measure of the asymmetry in the dot-lead couplings 
and Tr. Bias and magnetic field are converted to units of effective temperature, where 
T^ff = eV + Tb ase and Tf$f = Igl/^B-B + T&ase, allowing for a quantitative comparison between 
the response to the different perturbations. (The standard value of 2.0 was used in the 
conversion to units of T^^.) 

Adding the parameter Gb to Eqns. ([5]) and ([6]) allows for fits to data of G(T), G(V), 
and G(B). The confluence of these three relations leads to a value of Gb that most nearly 
produces a single energy scale, TV, for the three fits corresponding to a particular device. 
Thus the optimum Gb produces the best simultaneous fits for a single Tk- Certainly the 
parameters Go and Tk will depend upon the magnitude of the parameter Gb'-, thus, the value 
of Gb will also impact the values of the transport coefficients. In Fig. |2b G(T), G(V), and 
G(B) from sample A are plotted, along with their respective phenomenological curves, using 
T K = 7.32A" and G b = A.7e 2 /h. 

Analysis of 9 samples exhibiting a single channel spin-^ Kondo effect resulted in the 
following average values: ct = 4.97 ± 0.9; cy = 0.74 ± 0.14; cb = 0.89 ± 0.16 found via 
fits to Eqns. (SI)- (ED for energies up to T /Tk < 0.1. These values may be compared to the 
theoretical predictions of ct = 5.38, cy = 0.82, and cb = 0.55, in accord with the definition 
of Tk- We also find that the effect of magnetic field is well described by the theoretically 
predicted curve for T/T K < 0.1. We stress that Eqns. dlj-® are strictly valid only for the 
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FIG. 2: (color online), (a) Conductance in sample A vs. bias for several value of external magnetic 
field. Black line represents a fit to data using an equation consisting of a sum of two Fano peaks 
plus a cubic background (see Suppl. Info). Inset: Traces of dl/dVvs.V acquired at many values 
of B. Colorbar indicates magnitude of dl/dV ranging from 4.8e 2 //i (dark blue) to 5.0e 2 //i (dark 
red), (b) Equilibrium conductance as a function of temperature (green squares); zero-bias, base 

-eff 

same Gf, is subtracted from all data shown. Gray line is a fit to G(T) data using Eqn. Q. Black 
dash-dot line is a fit to G(B) data using Eqn. (0). Black line is a fit to G(V) data using Eqn. ([5]). 
The same value of Tk is used for all three fits. Horizontal gray dash-dot line signifies \Gq. 



temperature conductance as a function of effective temperature Tf^ (blue diamonds); and zero- 
field, base temperature conductance as a function of effective temperature TYtt (red circles). The 



particle- hole symmetric spin-^ Kondo model. To check for the internal consistency of our 



analysis, we compare with the method of Ref. 15] generalized to include finite magnetic 
field. From fitting the data of a representative sample to Eqns.(l)-(3), we find or = 4.91, 
cy = 0.72, c B = 1-49 and G(0, 0, 0)/G = 0.35. The large discrepancy between the two values 
for cb suggest that the g-factor of the device deviates from the assumed value g = 2.0. The 
renormalized parameters of the underlying Anderson model, as denned in Ref. 15[, turned 
out to be u = 0.68, q = 0.95 and £ = 0.263. While the value q = 0.95 points to a 
large degree of particle-hole asymmetry, we note that according to Fig.2a of Ref. [lsj], ct at 
u = 0.68 is hardly affected by it. 

Derivations of the finite bias conductance behavior at non-zero temperature have been 



worked out in the context of both the partic 



model 26M29| and the Kondo model 14 
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lole symmetric Anderson single level impurity 



321 ] . A scaling function for the conductance out 
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FIG. 3: (color online), (a) Conductance as a function of V for sample B at T/Tk < .65. (b) 
Scaled conductance, AG/ay, vs. (eV/kBTx) 2 for the same device at several temperatures, where 
AG = (1 - G(T,V)/G(T,0)), ay = c T a/{\ + c T {l~ l) (T/T K ) 2 ). We used a = .11 and 7 = 
.67. The dashed black line represents the associated universal curve. The data conforms well to 
the scaling function below (eV/fcB^k) 2 ^ -45 for all temperatures plotted. (c),(d) Normalized 
conductance plotted for 4 different samples. G(V) and G(B) data are plotted separately for 
clarity. Equilibrium conductance as a function of temperature (green squares) and (c) zero-field, 
base temperature conductance as a function of effective temperature T^ff ( re< ^ chcles), or (d) zero- 
bias, base temperature conductance as a function of effective temperature T^ f (blue diamonds). 



Same y-axis for (c) and (d). 



of equilibrium was outlined in Ref 16[ based on a low energy expansion of a form suggested 



by Ref 



22|, 



G(T, V) = G(T, 0) 



CtO. 



\ 



— (— V 



(7) 



1 + OT&-1) 

where the coefficients a and 7 represent the zero-bias curvature of the Kondo resonance 
peak as a function of V and the rate at which the resonance peak broadens and decays in 
amplitude as a function of T, respectively. Using a value for Tk inferred from Eqn. (pD), some 
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groups probing the Kondo state in nanotubes and quantum dots estimated 0.1 < a < 0.25 



23j | . Analysis of our new samples find a = 0.13 ±0.018 and 



and 0.5 < 7 < 1.65, Q, Q, 
7 = 0.89 ± 0.22. Note, however, that 7 depends on the lead-to-dot coupling asymmetry £ 
even in the particle-hole symmetric spin-| Kondo model. 15] 

In a past experiment measuring conductance in Ti/Au break junctions using two types of 
molecules two of the authors found values of a ~ 0.051 and 7 ~ 0.107. 17j A number of phys- 
ical influences were postulated as potentially affecting the devices under study leading to the 
exceptionally small coefficient values. For instance electron-phonon coupling could lead to a 



renormalization of estimated parameters. 
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34| Indeed, more recent theoretical studies have 



explored the impact of non-ideal conditions like finite on-site repulsion, tunneling barrier (i.e. 
molecule-lead coupling) asymmetry, particle-hole asymmetry, and charge fluctuations related 



to the possible proximity to the mixed-valence regime. 
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2n 
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35j They generally conclude 



that these conditions may impact a by lowering its value in the range of ~ 0.1 — 0.15. In the 



case of Ref . 



291 ] the authors anticipate 0.075 < a < 0.3. Employing a limited version of the 



treatment described here (field response is not available for the previously recorded data) to 



the data of Ref. ([17J) in order to account for background conductance, leads to the following 
coefficient values: a = 0.104 ± 0.01; 7 = 0.703 ± 0.12; c T = 5.25 ± 1.0; c v = 0.69 ± 0.18, 
which are in closer agreement with previously published results. Note that while ct is evi- 
dently fixed by the choice of s = 0.21 in Eqn. (j3J) and it is anticipated that cy = acx, the 
values found in our experiments differ slightly. This is due in part to the fact that Tk is not 
established by Eqn. (jlj) alone, and in part because of the more restrictive fitting range for 
the c's. For a and 7 the fits were limited to T/7V, e^/ZcsT^-, \g\fisB /ksTx < 0.2 to 0.3, 
depending on the sample. 

We have measured conductance in break junction devices in the low energy equilibrium 
and non-equilibrium Kondo regime. In some cases the peak conductance exceeds the unitary 
limit, highlighting the need for parallel transport channels to be taken into consideration. 
The use of multiple analytical expressions for the dependence on temperature, bias, and 
magnetic field provides a more reliable means of extracting the emergent energy scale, Tk, 
and it enables us to converge upon the proper conductance offset, which allows for the sepa- 
ration of background conductance from the cotunneling exchange process associated with the 
Kondo effect. We have applied this method to a group of new samples and found transport 
coefficients consistent with previous experiments and with theoretical expectations, partic- 
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ularly when non-ideal conditions are considered. When a limited version of this technique 
is applied to older data it also brings coefficient values in-line with expectations. Further 
inconsistencies may be addressed with additional considerations, for example accounting for 
a background conduction with a weak temperature or bias dependence. 

GDS would like to thank R. Willett and M. Peabody for technical support, and T.-C. Hu 
and R. Kopf for assistance with sample fabrication. DN acknowledges NSF DMR-0855607. 
SK and EM would like to thank T. Costi, and acknowledge support by the Comision Nacional 
de Investigation Cientifica y Tecnologica (CONICYT), grant No. 11100064 and the German 
Academic Exchange Service (DAAD) under grant No. 52636698. 
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SUPPLEMENTARY INFORMATION 



I. NON-SPIN-i KONDO EFFECT 



24 



thus consid- 



It is possible to find Ceo with spin S — 1 in break junction devices, 
eration was taken in this experiment to identify the spin state of the impurity in order to 
ensure that the analysis was applied only to samples containing a magnetic moment with 
S = \. The equilibrium temperature dependence of the Kondo resonance described by 
Eqn. (4) in the main text pertains only to the case of a spin-| impurity coupled to a single 
conduction channel. More generally, the following phenomenological function, connecting 
the low (T <C Tk) and intermediate energy (T ~ Tjc) regimes, can be used to describe the 
characteristic response of the system with respect to temperature, 



G{T) = G(0) 



-a s 



(8) 



Here T' K = T^-/(2 1/ ' Qs — l) 1 /^ 5 while as and are constants that depend upon the spin 
of the magnetic impurity and the number of conduction channels that may participate in 
screening the impurity. For an impurity with spin state | coupled to one channel we use 
a s = 0.21 and £ s = 2. 

To distinguish between a spin-| fully screened Kondo effect and a spin-1 underscreened 
(or fully screened) Kondo effect we attempted to fit the respective versions of the empirical 
form to the equilibrium data. A system is said to be fully screened when 2S = n, where S is 
the spin state and n is the number of screening channels. If 2S < n then the system is said 
to be underscreened. Because our devices effectively possess only one screening channel, we 



were primarily interested in the possibility of an underscreened effect. 37| The G(T) data 



was fit to Eqn. ([8]) using as = 0.33 and £5 = 0.89. These values were found via fits to the 



numerical renormalization group analysis in Ref . [36j . In most cases this produced very poor 
fits compared to the equation for the spin-| impurity, as shown in Fig. le of the main text. 
For the sample in Fig. SH]the fit of G(T) for S = \ and S = 1 produced roughly equal quality 
fits. Additionally, the magnetic field dependence of the base temperature conductance did 
not follow the expected format associated with a spin-| impurity. Due to the ambiguous 
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FIG. 4: (a) G(T, 0,0) vs. T for a sample that does not exhibit clear spin-^ characteristics. Gray 
line is a fit to the empirical expression in Eqn. (2) of the main text. The black dashed line is a 
fit to the analogous equation for an underscreened spin-1 Kondo effect shown in Eqn. (SI). The 
quality of the fits is similar, (b) G(0, 0, B) vs. B for the same sample. Of particular note is the 
increase in zero bias conductance as the magnetic field is increased up to approximately 2.6 K in 
units of effective temperature. While we were not able to resolve two separate peaks in traces of 
dl/dV vs. V acquired on the left side of the vertical gray dashed line, this behavior is consistent 
with a spin-1 Kondo effect in which the magnetic field drives a transition between a singlet and 
a triplet ground state. Inset: Energy as a function of magnetic field. Vertical gray dashed line 
corresponds to a magnetic field value that results in a degeneracy between the singlet |0 0) state 
and the triplet |1 —1) state. 



nature of this sample, and two others with similar properties, it was not included in the 
scaling analysis discussed in the text. 



II. DATA FITTING 

Equations 4 through 7 in the main text describing the conductance response 
to different perturbations, are valid only up to at most intermediate energies 
(T/Tk, eV/ksTK, \g\fisB /IzbTk ~ 1). As the perturbations are further increased relative to 
Tjc, it is expected that the measured conductance will deviate from these expressions due to 
the onset of transport processes unrelated to Kondo-assisted tunneling, yet this transition is 
not expected to show any clear demarcation. In the case of finite magnetic field (Eqn. (7)), 
it is expected that the conductance data will exhibit a stronger deviation from the analytical 
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formula since the theory is perturbative in the field, as opposed to non-perturbative results 
derived from numerical renormalization group techniques. 

To establish a consistent means of determining the low energy zero bias conductance 
for all traces of dl /dVvs.V a parabola is fit to the central portion of the resonance. It's 
maximum value is taken as the zero bias conductance. For traces in which the resonance 
peak is found to be offset from V = mV, commonly by < ±0.5 mV due to amplifier drift, 
the peak is re-centered for the ensuing analysis. 

The uncertainty in the coefficient values we found, as well as their corresponding values 
of Go and Tk, arise from a few different sources of error. Measurements of dl/dV at each V 
data point included some random noise, which varied in magnitude from device to device. 
There is some uncertainty in the thermometry calibration and temperature stability during 
the course of acquiring a trace of dl/dVvs.V. Finally, the range of data points included 
in fits can impact the value of the fitting parameters. This is especially noteworthy for 
conductance probed as a function of temperature. The comparatively limited control and 
lower number of data points collected as a function of temperature, in contrast to what was 
acquired with respect to bias and magnetic field, makes the extraction of Tk susceptible to 
significant differences for small changes in the fitting range. 

To maintain the validity of our analysis the fit ranges were constrained to lower energies 
so as to avoid fitting to data outside the applicable range of the functions. For a given 
sample the data range used for the fits was approximately equal for G(T) and G(V), and 
was commonly somewhat smaller for G(B) due to the limited range of validity of Eqn. (6), 
as mentioned above. For a given sample small differences between the range of G(T) and 
G{V) occur because of the larger spacing between data points for the measurements of G(T). 
The quality of the fits was compared on the basis of the root mean square deviation of the 
data from the analytical expressions. For G(V) and G(T), fits were performed up to T /Tk 
and T^j/Tk < 0.2 to 0.4 depending on the sample, and for G(B) fits were performed up to 
T^jj/Tk < 0.08 to 0.15 where the numbers were also sample dependent. 

In accord with the Zeeman effect, the application of an external magnetic field, B, lifts 
the spin degeneracy of the many-body Kondo singlet, and results in the splitting of the 
zero-bias resonance peak. As B is increased from T to some field value Be, the Kondo 
resonance, dl/dV vs. V, is observed to broaden and decrease in height. For B > Be the 
splitting is large enough that two peaks can be resolved. Their position and amplitude is 
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ascertained from by a fit to the sum of two fano lineshapes plus a cubic background, |3[ 

-B\V\\ (9) 



dl A (ei + gi) 2 , . (62 + g 2 ) 2 

— = y + Ai— — — + A 2 —— — = 
dV 1 + ei 2 1 + e 2 2 



where €1,2 = (V — xx^/wx^, £1,2 is the center of a peak, and w\p, is its half- width at half 
maximum. This is demonstrated in Fig. 2a in the main text, in which the solid black line 
is a fit of Eqn. (SEJ) to the dl/dV vs. V data corresponding to B = 8T. 

III. ZERO BIAS ANOMALY STATISTICS 

Thermally-assisted electromigration was performed on 62 break junction samples. Of 
the 35 samples decorated with C6o> 12 exhibited some sort of zero bias anomaly, but only 
9 displayed clear signatures of the single channel spin-^ Kondo effect. Molecules were not 
deposited on 27 samples. These were utilized as control devices in order to test for the 
likelihood of observing zero bias anomalies unrelated to transport through a Cgo molecule. 
In this control group a zero bias conductance peak was observed in two devices. The rate 
with which we expect to see Kondo-related transport are commensurate with the findings 
of past investigations. 0, Isl, tl7 [ 

The rate of occurrence of a zero bias conductance peak in samples decorated with Cgo 
is much greater than that of the control devices. This leads us to believe that transport in 
these samples is related to the Ceo- However, we acknowledge the possibility that transport 
leading to the resonance peak in the samples identified as S = | Kondo effect due to a 
single unpaired spin situated on a Ceo molecules could be due to an alternate process or 
composition. 

Another experimental realization of the Kondo effect has been demonstrated in chemically 



homogenous magnetic nanocontacts without any clearly defined impurity. In Ref. [12] the 
authors argue that a small portion of one contact may function as the local moment. This 
local moment could form in the contact atoms because of their smaller coordination and 
reduced symmetry, compared with the bulk. The metal used in our break junctions is 
Pd. While Pd is paramagnetic in its bulk state, only a small change in band structure is 
required to trigger itinerant ferromagnetism. Nanoscale structuring via the electromigration 
technique will affect grain boundaries and lead to undercoordinated atoms at the break 
junction site, which can alter the density of states at the Fermi level thereby affecting the 



14 



local magnetization. 38N40 1 



Our study is primarily focused on probing the Kondo state, regardless of what constitutes 
the local impurity. Therefore, the possibility that transport may be occurring through a local 
moment other than a single Cgo molecule does not negate the consistency of the results. 
Furthermore, part of this study is aimed at contrasting the results obtained from different 
experimental realizations of the Kondo effect, and a magnetic impurity in a break junction 
device, regardless of its makeup, constitutes a fundamentally different system compared to 
experiments probing Kondo scattering in semiconductor and nanotube structures. 
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